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Abstract 

Using the relation between the harmonic oscillator wave function and the light 
quark scattering form factor, the expectation values of colour-straight four- 
quark operators are evaluated and found to be directly proportional to the 
cubic power of the oscillator strength. It is predicted that the ratio r(Ab) /t(B) 
~ 0.79(0.84) due to the factorizable (nonfactorisable) piece, against the exper- 
imental 0.79 ± 0.06. Notwithstanding the numerical prediction, the present 
study shows that the four-quark operators play a role as far the lifetimes of 
b-flavoured hadrons are concerned. 
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I. INTRODUCTION 



In the description of heavy hadron decays by heavy quark expansion (HQE), the 
preasymptotic effects appearing at next-to-leading order and beyond are vital in predict- 
ing the decay properties accurately. These effects are due to the operators of dimensions D 
= 5 and 6. At I/itlq, the operators are suppressing the leading order. The evaluation of 
D = 5 operators which describe the motion of the heavy quark inside the hadron and the 
chromomagnetic interaction is definite. The estimation of the D = 6 operators which are 
four-quark operators (FQO) containing both heavy and light fields is based on the vacuum 
insertion assumption for mesons and on the quark models for baryons. Though their effects 
are negligible as the heavy quark's volume occupation is of the order of (Aqc^d/ttiq) 3 but 
are considerably enhanced due to partial compensation by the four-quark phase-space, these 
operators are predicting the lifetime differences in the world of the hadrons of given flavour 
quantum number. Therefore, accurate value of the FQO is necessary due to the confronta- 
tion existing between theory and experiment over the hadronic properties like the exper- 
imentally smaller than theoretically expected lifetime of A& and the theoretically smaller 
than experimentally predicted semileptonic branching ratio of B-meson. Theoretically upto 



order (l/m fe 2 ) §g 




1 + 0.01 



t(B°) 
r(A>) 
t(B°) 



0.9 



(1) 



whereas experimentally ||, 



r(B+) 



1.04 ±0.04 



t(B°) 
r{B s ) 



0.99 ±0.05 



r{B°) 



0.79 ±0.06 



(2) 
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The agreement among the B mesons is as expected of but not so between the A& and 
B. The later issue continues to be central point of the physics of heavy quark hadrons. It 
is suspected that the explanation for these discrepancies within the HQET framework is 
hidden in the yet-not-satisfactorily-understood FQO. 

As regards the evaluation of the FQO, there were two works which attempted 
to explain substantially the enhanced decay rate of A&, whereas the work of P. Colangelo 
and F. De Fazio f| which is QCD sum rules based prediction leads to conclude that the 
reason for the smaller lifetime is not due to FQO. In Ref. Q, the authors evaluated the 
FQO parameterising the matrix elements in terms of hadronic parameters which are not 
practically known. But the parameters have been calculated using QCD sum rules 
However, the prediction is not able to account for the lifetime difference between A& baryon 
and B meson. On the other hand, the author of Ref. |J used quark model and accounted for 
the FQO for 13% of the required enhancement in the A& decay rate. The above estimation 
used yet to be confirmed result of DELPHI collaboration || on the mass splitting of and 
Eft. The same method has been modified by taking the logarithmic dependence of the wave 
function at the origin and this explains the difference between the decay rates of B meson 
and A b baryon to the extent of 40% 0. Since the striking point in the evaluation of the 
FQO is not yet obtained to clear the situation in one way or the other, it is important and 
interesting to explore other avenues to estimate the four-quark matrix elements. 



In this paper, we adopt the colour-straight formalism approach of [[L0) to evaluate the 
expectation values of the four-quark matrix elements. On the specific choice of the harmonic 
oscillator wavef unction model for the form factor and slightly different potentialr for meson 
and baryon, it is found that 

r(B+) 

1.00(1.03) 
1.00(1.02) 



t(B°) 
r{B s ) 



t(B°) 

- = 0.79(0.84) (3) 



r(A b ) 



t(B°) 
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where the values given within brackets are due to non-factorisable part of the FQO. These 
values are in agreement with the data, eq. (2). 

In a recent work flOfl , Pirjol and Uraltsev discussed the four-fermion operators on certain 
quantum mechanical basis. In the nonrelativistic quark theory, the wave function density 
and diquark density are related to the associated operator 

(knb^qjT^), (4) 

where the I\,j are arbitrary Dirac structures, through 

' < B\(bb)(qq)\B >= \*(0)\ 2 (5) 



2M B 



^- < A b \(bb)m\A b >= J d 3 y|*(0,2/)l 2 (6) 

for meson and baryon respectively The operators in Eq. (4) are colour singlet. As colour 
flows freely for these operators, they are called colour-straight. The expectation values of 
the operators in Eq. (4) are related to the observable, the transition amplitude of the light 
quark scattering off the heavy quark. Thus the determination is based on the knowledge of 
the light quark scattering form factor. 

The wave function at the origin, in momentum representation, is given by 

The transition amplitude is then the Fourier transform of the light quark density distribution: 

F(q) = 21^ < H ^M°)\ H b(°) >= J rf 3 xvT(x)vI/(x)e»> x (8) 
Integrating over all q yields the expectation value: 

/ (2^ F(q) = im{2 = 2M^ < H ^(0)\H h > (9) 

For any Dirac structure T, the light quark current density and the light quark transition 
amplitude are given by: 



J r (x)=gTg(x); A r (q) = JL < H b (q)\J r (0)\H b (0) > (10) 

where the Jr(0) is gauge invarinat operator and not required to be a bilinear. Thus, for 
spin-singlet operators, we have 



< H(q)\bbJ r (0)\H b (0) >= J — < H b (q)\J r (0)\H b (0) > (11) 
And, for spin-triplet operators, similarly 



< H(q)\ba k bJ T (0)\H b (0) >= J < S k H b (q)\J r (0)\H b (0) > (12) 

with S/2 being the b-quark spin operator and 

\S k H b >= J rf 3 x6a fc 6(x)| H b > (13) 

Equations (11) and (12) represent the general structure of four-quark operators. The above 
operators are local as required for by the HQE in the sense that the light quark operators 
enter at the same point as the heavy b-quark operators. These relations, Eqs. (11) and 
(12), hold equally for different initial and final hadrons having different momenta smaller 
than m&. 

Equation (4) resolves into spin-singlet and spin-triplet operators for 1^ = 1 and T b = 
775 ( = a ) respectively. The light quark elastic scattering is described by the form factor 
F(q 2 ). In Ref. [0, the exponential ansatz and the two pole anstaz are used for the form 
factor. Both of them lead to a determination which make no difference for meson and 
baryon. This cannot truly be the case to represent ^(O)) 2 , the measure of the expectation 
values of FQO, for a baryon and a meson. 

In the next section, the choice of the representation, harmonic oscillator wave function 
for the form factor and the potentials are discussed. The evaluation of the expectation values 
of the factorisable part of the FQO is presented in Sec. III. Corresponding nonOfactorisable 
part is given in Sec. IV. Estimation of the lifetimes ratio and the conclusion are given in 
Sec. V and VI respectively. 
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II. HARMONIC OSCILLATOR WAVE FUNCTION MODEL FOR FORM 

FACTOR 



As will be discussed in the following sections, the expectation values of the colour-straight 
operators are parameterised in terms of a single form factor for both B-meson and baryon. 
The extraction of the form factor involves assumption of a function such that it satisfies the 
constraints on the form factor that F(q 2 = 0) is equal to the corresponding charge of the 
hadron. Then the form factor has to be extrapolated into the region where q 2 > 0. We take 
the hadronic wave function of ISGW harmonic oscillator model [111] f° r the f° rm factor. 

The wave functions of ISGW model are the eigenfunctions of orbital angular momentum 
L = satisfying the overlapping integral 

/(q) = J r 2 dr^/* f (r)^/i(r)j (qr). (14) 

The overlapping integral can be equated to the form factor. Hence for diffrent initial and 
final hadrons 

F( q 2 ) = N 2 exp[-q 2 /2(Pj + Pf)] (15) 

where N is the normalisation constant given by [2f3f(3i/((3j + /3 2 )] 3 / 2 and /3's are oscillator 
strengths. For same initial and final hadrons, the transition amplitude is 
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= ^ = l*W (16) 

From the above equation, which is the central point of dicussion of this paper, it is obvious 
that the transition amplitude and hence the expectation values of four fermion operators 
are proportional to the third power of the oscillator strength of the hadron. 

The calculation of /3's can be made using the QCD inspired potential. In the present 
calculation, we use the potential for B-meson containing the Coulomb, confining and a 
constant terms: 

V(r) = ^- + br + c (17) 



For a = -0.508, b = 0.182 GeV 2 and c = -0.764 GeV, and for quark masses m g = 0.3 GeV 
(treating m„ = m d ), m s = 0.5 GeV and m b = 4.8 GeV, using variational procedure with the 
wave function given in position space as, 

= ^ r '< 2 («) 

(3 Bq = 0.4 GeV and (3 Bs = 0.44 GeV. 

For A;, baryon, following the similar procedure but for the potential of the form 

V(r) = ^(^- + br + Pr 2 + c^ (19) 

where the r 2 term is a harmonic oscillator term justifying the consideration that A& be a 
two body system and the same wave function of Eq. (18), one gets (3^ b = 0.72 GeV and for 
Ps b = 0.76 GeV for the values of mass of the diquark system 0.6 and 0.8 GeV respectively. 
The large value for the /3\Lambda is due to the presence of the 0(r 2 ) term in the potential. 
Otherwise, the value of (3a is no different than that of B s . These values are used in the 
subsequent calculations in this paper. A comment is in order on the choice of the same wave 
function for baryon as for meson: In the usual procedure, the ground state wave function 
for baryon is 

Aground « e^*^ 2 (2Q) 

where r A p are the internal coordinates for three body system. Due to the idea of considering 
the A b as a system containing the bound state of light quarks and a heavy quark, the 
separation of the two light quarks which make the bound state is treated negligibly. This 
allows then that the baryon is a system of two body. It is a reasonable approximation only. 
The difference between a meson and a baryon is essentially due to the value of the oscillator 
strength. 

III. EXPECTATION VALUES OF THE COLOUR-STRAIGHT OPERATORS 



We evaluate the expectation values of the colour-straight operators only for the vector 
and axial-vector currents. Nevertheless the other currents can also be studied in the same 

7 



fashion. Both the currents are possible for B-meson while axial currents vanish for A& baryon 
due to the light degrees of freedom which constitute a spinless bound state. 

Essentially there is no difference between the exponential ansatz and the harmonic os- 
cillator wave function in representing the behaviour of the form factor but they differ while 
fixing the scale: in the former case, the hadronic scale of one GeV is used whereas in the 
latter the same has been fixed by solving the Schrodinger equation. The two pole anstaz is 
based on the well founded experimental values. Basically the use of the harmonic oscillator 
wave function of the constiuent quark model is an alternate picture but in the very same 
lines of the two ansatz. 

Hereinafter the operators are referred to by the following notaion: for meson 



< O q VA > H =< H\bT VA bqT VA q\H >; < T VA > H =< H\bT VA t a bqT VA t a q\H > (21) 



where T v ,a = 7^, 7/i75 and t a t b = 1/2 — l/2N c . In what follows, q stands for u and d quarks 
and s for strange quark. And < Oy,A >, < Ty,A > will be respectively denoted as u>v,a, t v,a- 
For baryon, < 0(T) V > correspond to Ai j2 respectively. 



The parametrisation of the matrix element of the colour-straight operators for vector 
current is 



with the constraints F B (0) = 1 for valence quark current and F B (0) = for sea quark current. 
The former is relevant for the b-meson composition of quarks bq. Then the corresponding 
transition amplitude is 



A. B-meson 



< B(q)\J Tv \B(0) >=- VfM F B (q 2 ) 



(22) 




< B(q)\O 9 v \B(0) >=-v, 



(23) 



Under isospin SU(2) symmetry, 
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< B~\Oy\B~ >= -2.88 x 1(T 3 GeV 3 ; < B~\T$\B~ >= -9.61 x lO~ 4 GeV 3 (24) 
For B s , we have 

< B s \O q v \B s >= -3.83 x 1(T 3 GeV 3 ; < B S \T V \B S >= -1.28 x 1(T 3 GeV 3 (25) 

If the case of violation of isospin symmetry, and SU(3) flavour symmetry, is considered, then 
there comes Cabibbo mixing of eigenstates for d and s quarks among themselves. That is, 
d! = dcos9 c + ssinO c and s' = scos8 c + dsin9 c . This we do not consider here. 

For axial current there are two form factors which are related to one another due to 
conservation of the axial current, d^J^ = 0, in the chiral limit. By the Goldberger-Treiman 
relation |T2| which equates axial charge form factor to the coupling gB*Bn at q 2 = 0, the 
operators involving axial-currents are estimated in terms a single form factor. Thus, given 
the value of the coupling g, the extraction of the transition amplitude is similar to the 
B-meson case. 

Making use of 

(S b e)\B(q)>=\B*(q,e)> (26) 
and Eq.(12), the expectation values for the axial vector currents are given by 

< B(q,ey\ £ q W q\B(0) >= e^GfV) - (e* q)q,G ( §\q 2 ) (27) 

q=u,d,s 

< B(q,ey\q\ a llll5 q\B(0) >= {e^Cf V) - (e^GJV)}^,- (28) 
Finally the following equality leads the absence of the structure (e*g)v M 

< B*{q, 6)1^(0)15(0) >*=< B(0)\j, 5 (0)\B*(q,e) >=< B*(0, e)|^(0)|B(g) > (29) 
Following the Goldberger-Treiman relation, we have 

G 1 (q 2 )=q 2 G (q 2 )=ge^^ 2 (30) 
Correspondingly, the expectation values ars 
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< B-\0\\B- >= 8.63 x 1(T 5 Gel/ 3 ; < B~\T A \B~ >= 2.88 x 1(T 5 GeV 3 (31) 
< B s \O s A \B s >= 1.15 x 10" 4 Gel/ 3 ; < fi s |T||fi s >= 3.84 x 10~ 5 GeV 3 (32) 



We have taken in the above esitmates the value g = - 0.03 [13,14 . 



B. A b baryon 

For A b baryon, treating u and d quarks equally, 
< A 6 |0^|A fe >= -1.69 x 10~ 2 Gel/ 3 ; < A b \T$\A b >= -5.64 x 10~ 3 GeV 3 (33) 
In the case of E b , we have, 

< E b \ °v\-b >= -2.01 x 10~ 2 GeV 3 ; < E b \ ]T O q y\E b >= -6.72 x 10" 3 GeV 3 

q'=u,d,s q'=u,d,s 

(34) 

There are corrections additionally to form factors due to charge radius. The same can be 
ignored as we are looking at the wave function density at the origin. 

IV. NON-FACTORISABLE PART OF THE FQO 

The nonfactorisable part of the FQO come in four. The following is one of the ways of 
parameterising them [[|. 

1 



2M B 
1 



< B\{bq) v _ A {qb) V -A\B >= f B M B B x /2 (35) 



2M B 
1 



2M V 
1 



< B\(bt a q) v _ A (qt a b) v _ A \B >= f 2 B M B e 1 /2 (36) 



< B\{bq) s „ P {qb)s-p\B >= f B M B B 2 /2 (37) 



< B\{bt a q) s .p{qt a b) s _ P \B >= f B M B B 2 /2 (38) 



2M B 
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where -81,2 and ei_2 are hadronic parameters. They are related to wy,A and tv,a which are 
the expectation values of the operators Oy,A and T Vt A as defined earlier. 

P B M B B l = (t> 1 = A(Tv + T A ) + — (u v + u A ) (39) 
f B M B B 2 = 2 = -2(t v - t a ) - -^-(tov - u A ) (40) 



flM B e x = Pl = -u-(t v + ta) + (1 - T7j)K + ua) (41) 

f B M B €2 = p 2 = -^(r v - T A ) - 1(1 - ^2 )K - wa) (42) 
In the case the A& baryon, the nonfactorisable piece corresponds to 

< (&9)v-a(?&)v-a) = -7TTT A i - A 2 (43) 

Zl\c 



V. DECAY RATES AND LIFETIMES 



The decay rates of the b-flavoured hadrons are given by 

Ai + 3A 2 



T(H b -> HM) 
where, with x = m 2 c /m 2 



G)ml\V cb \ 2 r 

1927T 3 



1 + 



2m? 



)/(*) - 



2mn 



(44) 



/(x) = 1 - 8x + 8a; 3 - x 4 + 12x 2 lnx 



(45) 



f(x) = (l-a;) 4 (46) 

are the QCD phase space factors and Ai and A2 correspond to the motion of the heavy 
quark inside the hadron and the chromomagnetic interaction respectively. These values are 
taken to be -0.5 GeV for meson and -0.43 GeV for baryon and 0.12 GeV for meson. The 
chromomagnetic energy is zero for all baryons except Qq. 
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Equation (44) is further supplemented by the FQO at the order (l/nig) in the HQE. At 
this order due to the light quarks there are three processes: Pauli interference (PI), Weak 
Annihilation (WI) and Weak Scattering (WS). The PI plays a predictive role in the charged 
meson and the A fe baryon. The PI becomes constructive at the tree level whereas it becomes 
destructive if radiative corrections are considered. The WS takes place in the neutral meson 
as well as A&. 



1. Lifetime ratio of B and B& 

Although the difference between the lifetimes of the charged and the neutral B-mesons 
is almost a settled issue, we check the once again using the expectation values of the colour- 
straight operators. This difference is attributed to PI and WA. Neglecting the WA as it is 
strongly CKM suppressed the result for the PI is 

Ar'(iT) = r 24 7 r 2 C <Q ^ >B ~ ~ 3 < °^ >B ~ (47) 

where C = c\ — c 2 _ + jf{c\ + c 2 _) and r = G/ ^ 2 ^ • ^ ne vames f° r Wilson coefficients 
are: c + = 0.84 and c_ = -1.42 with N c = 3. For r , m fe = 4.8 GeV and \V cb \ = 0.04 are used. 
Then the ratio is 

T(B = 1 00 (48) 

This agrees well with the one obtained in terms B-meson decay constants. 
The decay rates due to spectator quark(s) processes are: For B , 

Ar-f(B-) = ^\V cb \\l - xf[{2c\ - c 2 _)^ + 3(4 + cl) Pl ] (49) 

Hence the ratio is 1.03. 



2. Lifetime ratio of B and B s 



The lifetimes difference between the two neutral mesons B s and B^ is due to W exchange. 
The numerical result is 
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4^4 = 1-00 (50) 
r(B d ) 

Corresponding to the nonfactorisable part, we get the decay rate 

AP*(B.) = ^\V cb \ 2 (l - ^[ {2C+ - C - )2 ((1 - x)4>\ - (1 + 2xW 2 ) + (C++ 2 C - )2 ((1 - x)p\ - (1 + 2, 

(51) 

Therefore the ratio becomes 1.02. 



3. Lifetime ratio of A b and B 

In the HQE, the difference in lifetimes between mesons and baryons begins to appear at 
order l/nig. Nevertheless it is dominant at third power in 1/mg. At this order, the FQO 
receives corrections due to WS and PI. They are 



T ws (A b ) = 927r 2 r c 2 <Q ^ 3 >Ab (52) 

ml 

V PI {K b ) = -487r 2 r Ci <Q V >A '' ( 53 ) 

where C\ = — c + (2c_ — c + ). As mentioned earlier PI is destructive for radiative corrections 
and enhances the decay rate leading to smaller lifetime for A 6 . The effect of WS, on the 
other hand, is colour enhanced and its consequence is smaller. Hence, 

$H-™ <«> 

The deacy rate modified by the nonfactorisable piece is given by 

G 2 m 2 - 

Ar(A 6 ) = -^A[4(l - x)\cl - 4) - (1 - xf{\ + x)(c_ - c + )(5c + - c_)] (55) 

107T 

where A stands for the term in Eq. (43). Correspondingly, the ratio is 

™= 0M (56) 

In mesonic cases, the nonfactorisable piece gives a little bit higher values. In particular, the 
ratio of the lifetimes of the baryon and meson is significantly larger. 
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VI. CONCLUSION 



In this paper, we have evaluated the FQO for beauty hadrons. Though the spectator 
effects are suppressed by powers of (Agc^/mg) 3 , in the HQE for inclusive decays, they 
cannot be neglected. We have expressed the four-quark operators in terms of light quark 
scattering form factor which are in turn related to the harmonic oscillator wave function. 
The use of the wave function model is to replace the exponential and two pole ansatz used 
W\. Basically both are same. The distinction arises only due to /3, the oscillator strength 



in 



of the model. Interestingly this simple alternative predicts the lifetimes ratio of A& and B 
closer to the experimental value. 

On the other hand, the nonfactorisable part does not have much effect in the case of 
mesons. But it keeps still away the ratio between B and A& away from the experimental 
value. As far the B-mesons are concerned, the present study once again affirms the existing 
predictions. In this case too, there are omissions like SU(2) and SU(3) symmetry breaking. 
They may play a role but too negligibly. 

Finally, we conclude that we have taken one, which is dominant, of the sources of the 
preasymptotic effects and shown that it predicts the lifetime of the A& close to the experi- 
mental figure. As we have not taken into account all possible corrections to the four-quark 
operators, the present prediction can be considered at least indicative in order to look into 
the four-quark as well as six-quark operators more seriously. However given the basis pro- 
vided in [jToH , the prediction has to be believed. Of course, this prediction can be checked 



by lattice studies. A refined analysis of b and c flavoured baryon lifetimes will be published 
elsewhere. 
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